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We study how the formalism of the Hierarchical Reference Theory (HRT) can be extended to 
inhomogeneous systems. HRT is a liquid state theory which implements the basic ideas of Wilson 
momentum shell renormalization group (RG) to microscopic Hamiltonians. In the case of homoge- 
neous systems, HRT provides accurate results even in the critical region, where it reproduces scaling 
and non-classical critical exponents. We applied the HRT to study wetting critical phenomena in a 
planar geometry. Our formalism avoids the explicit definition of effective surface Hamiltonians but 
leads, close to the wetting transition, to the same renormalization group equation already studied by 
RG techiques. However, HRT also provides information on the non universal quantities because it 
does not require any preliminary coarse graining procedure. A simple approximation to the infinite 
HRT set of equations is discussed. The HRT evolution equation for the surface free energy is numer- 
ically integrated in a semi-infinite three-dimensional Ising model and the complete wetting phase 
transition is analyzed. A renormalization of the adsorption critical amplitude and of the wetting 
parameter is observed. Our results are compared to available Monte Carlo simulations. 



I. INTRODUCTION 



Effective interfacial Hamiltonians are widely used to describe the large-scale fluctuations that occur near surface 
critical phenomena such as wetting |j| . Of course these models are not truly microscopic but are usually considered 
valid for length scales larger than some appropriate cut-off. The prevailing belief is that interfacial models may be 
derived from more microscopic approaches if the bulk degrees of freedom are integrated out. Needless to say this is an 
extremely difficult task and interfacial models still retain a partly phenomenological status. Therefore it is interesting 
to develop a genuine microscopic approach to study interfacial behaviors. Effective Hamiltonians are used to describe 
the critical behavior also in homogeneous systems because they allow a direct implementation of renormalization 
group (RG) ideas. For such bulk systems, the Hierarchical Reference Theory |3| provides a systematic way to 
derive an effective Hamiltonian from a given microscopic model Q . This theory, which implements the basic ideas of 
Wilson momentum space renormalization group for microscopic Hamiltonians, allows to derive an exact hierarchy of 
differential equations describing the evolution of the free energy and of the n-point correlation functions of the system 
when fluctuations on larger and larger length-scales are included. This hierarchy of differential equations can be 
closed, for example, by imposing an approximation for the pair correlation function, usually of the Ornstein-Zcrnikc 
(OZ) form. Already at this level of approximation, the HRT shows genuine non classical critical behavior. 

The extension of the HRT formalism to inhomogeneous systems would enable us to study the interfacial phase 
transition starting from the microscopic model, without the explicit introduction of effective interfacial Hamiltonians. 
Given two homogeneous phases, a and f3, in contact with a third passive phase (a wall), we study the particular case 
of complete wetting in planar geometry |]|, i.e. the phase transition which corresponds to the growth of an infinitely 
thick liquid layer close to the wall when bulk phase coexistence is approached at fixed temperature. The standard 
fluctuation theory of the wetting transition, beyond the mean field approximation, is described by using the capillary 
wave (CW) effective Hamiltonian in d dimensions 



H CW = I 



|(VZ(x)) 2 + W/G(x)) 



(1) 



where ^(x) is a collective coordinate which represents the distance of the a — (3 interface from the wall, W(l) is the 
effective interface potential which describes the effective interaction between the wall and the a — (3 interface and £ is 
the stiffness coefficient of the interface. For isotropic models £ can be identified with the surface tension, but it may 
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also depend on the curvature of the interface. The form of effective interfacial potential W(l) depends on the range of 
the microscopic interaction. In the case of short range potentials, for large I it is parametrized as pi ml IsL Ifll ml Hl| : 

W{1) ~ -Aexp(-ml/0 + -Bexp(-nZ/f) (2) 

where to, n are two dimensionless constants with n > to. The parameter £ can be identified with the bulk correlation 
length. It is important to note that the use of such an effective surface Hamiltonian can be justified only when 
bulk fluctuations are small i.e. far for bulk critical point Q. In three dimension, renormalization group studies on 
the Hamiltonian 8| ITol ] predict that the wetting critical behavior is non-universal depending on the value of the 
wetting parameter: 



k B T 
47r££ 



2 (3) 



where is the true bulk correlation length which governs the exponential decay of correlations in real space 0, H3 . 
At the critical wetting transition (i.e. wetting at coexistence) the critical exponents depend on this parameter [8j 
while at complete wetting (i.e. wetting approaching coexistence) the critical exponents are predicted to remain mean 
field like but critical amplitudes are ui dependent. For example the interface height I grows as: 

i~(l + |)log(A M ) (4) 

where A/i = ji — /io is the deviation of the chemical potential from the value at coexistence. The mean field approxi- 
mation to an effective 4> A hamiltonian provides a value u m f — via Eq.(4). 

The three-dimensional semi-infinite Ising model is one of the simplest microscopic models which may be simulated 
to test the RG predictions. The value of the wetting parameter lj is a function of the temperature and has been 
theoretically estimated 0,0]: to ~ 0.8 for T r > T > 0.6T C . Extensive Monte Carlo 0] simulation studies appear to 
be consistent with w ~ 0.3 at T S i m ~ 0.663T C jT^] a much smaller value than the predicted one. One possible solution 
of this discrepancy between the simulations and the RG results is the inadequacy of the CW Hamiltonian and the 
introduction of more general effective interfacial Hamiltonians [TtI fl8L IT^. l2fj| . 

The HRT theory, suitably generalized to deal with inhomogeneous systems, can be used to study the effects of 
thermal fluctuations beyond the mean field behavior without reference to the CW approach and so without introducing 
u> as an external parameter. We first obtain the HRT surface evolution equation for the Ising model which is then 
numerically solved in the case of complete wetting in three dimensions. The divergence of the adsorption T is studied 
when coexistence is approached and the critical amplitude is evaluated and compared to simulations. This work 
is organized as follows: in Section II we extend HRT approach to the case of inhomogeneous systems with planar 
geometry. We derive the evolution equation for the bulk and surface free energy. We also show that in the asymptotic 
region the HRT equations reduce to the known RG approach. In Section III we study complete wetting for a lattice 
gas model with nearest- neighbor interactions (which is equivalent to the Ising model). We investigate how fluctuations 
modify the mean field picture by integrating the HRT surface equation. In Section IV we briefly summarize the most 
relevant results and compare them to the available simulation data. 



II. HRT FOR INHOMOGENEOUS SYSTEMS 



A. Evolution equations 



The starting point in the derivation of HRT equations is the separation of the inter-atomic potential v(r) in two 
parts: 

v(r) = vji(r) + w(r) (5) 

where vn(r) is the short-range repulsive part of v(r). The thermodynamic and structural properties of the system 
with interaction vr(t), the "reference system", are considered known, at least numerically. It is also assumed that 
there is no phase transition in the reference system. Instead w(r) is a (mostly) attractive term, which triggers the 
liquid-vapor phase transition. Using this separation and performing a Legendre transformation on the grand partition 
function, a formal diagrammatic expansion for the Helmholtz free e nerg y can be written to all orders in perturbation 
theory 0] • We implement the basic ideas of Wilson's RG approach |21j within such a formal perturbative expansion 
to study how the bulk and surface thermodynamic quantities evolve due to the inclusion of fluctuations. This can 
be done by introducing a sequence of intermediate potentials characterized by an infrared cut-off in Fourier space, 



depending on a parameter Q. In the HRT for (off-lattice) bulk systems this cut-off is spherical in Fourier space, 
thereby respecting the isotropy of the interaction. Here we want to consider the case of inhomogeneous systems in 
the presence of a planar wall, where only cylindrical symmetry survives in wave vector space. Moreover, in wetting 
phenomena the correlation length parallel to the surface £|| and the correlation length perpendicular to the surface 
are related by Q, : 

f[. ate)] ^=3 

{(,• d<3 

This relationship suggests that the fluctuations perpendicular to the wall diverge much more slowly than the fluc- 
tuations parallel to the wall. Therefore it is natural to define a sharp cylindrical cut-off |'2'2t which prevents long 
wavelength CW critical fluctuations. Within the HRT approach, this is implemented by defining the sequence of 
intermediate potentials: 

-»(«,.) m 

where k is the component of the wave vector parallel to the surface and q is the component normal to the wall. The 
system characterized by the potential v®(r) — vr(t) +w®(r) will be named Q-system. For Q — > oo, reduces to vr 
while for Q — > the full interaction is recovered. Therefore it is natural to look for evolution equations governing the 
change in the physical properties of the Q-system as the cut-off is varied. In order to derive such evolution equations, 
we consider the structure of the perturbative series which defines the free energy of the model. In the perturbative 
diagrammatic expansion of the free energy of the Q-system A® , every loop contains one or more cf> — —flw® bonds . 
Therefore, due to the vanishing of w®(k, q) for k < Q, A® is defined by the same perturbative expansion as the full 
free energy of the model, A, where each (longitudinal) momentum integration is limited by an infrcred cut-off Q. The 
correspondence between a cut-off in the potential and a cut-off in fluctuations (i.e. in the momentum integrations) 
is valid with the single exception of the first, mean field, diagram of the perturbation series which does not contain 
any loop. This contribution is in fact discontinuous in Q because it contains just the n — Fourier component of the 
attractive potential 10 which is zero for every Q ^ and finite for Q = 0. However we can introduce a modified free 
energy A® which is continuous in Q and is simply related to A® 

- / 3^ = - /3 ^-I[0(r = O)-0 < 3(r = O)] / d d v p(r) + 1 /d d nd d r 2 [^(n, r 2 ) - ^(n, r 2 )] p(n)p(r 2 ) (8) 

Analogously, we can introduce the direct correlation functions , for each Q-system by functional derivation of the 
free energy A® with respect to the local density p(r). These correlation functions are continuous in Q for n > 3 
because there are no zero loop contributions in the corresponding perturbative expansion, but the two point direct 
correlation function is discontinuous at k — Q. So we introduce a modified function which is continuous, being the 
second functional derivative of A® with respect to the density: 

C 2 Q («,«) = C 2 Q (q, k) + [(f>(q, k) - <j> Q (q, «)] ■ (9) 

It is apparent from the definitions © and JjJJ that A® and coincide with the free energy and direct correlation 
function of the fully interacting system respectively when Q — > 0. On the other hand in the limit Q — > oo these 
modified quantities reproduce the mean field approximation for the free energy and the direct correlation function, 
contrary to A® and which reduce to the reference system quantities. This suggests that the HRT procedure does 
indeed describe the growth of fluctuations on top of the mean field approximation, as in the RG approach. 

The simple relationship between the cut-off-dependent modified quantities (.4^ and C§) and the Q-system properties 
(A® and C^) allows to derive the evolution equations describing how A® and change when Q is decreased from 
infinity to zero, i.e. when fluctuations on larger and larger length-scales are included. The perturbative expansion 
of the free energy can be specialized to the case where the reference system is the Q" s y s te m and the perturbation 
potential is 

Sv%,n) = v^{q,n) - v%, K ) = { ^ < K < Q (10) 



where SQ > is an infinitesimal shift in the cut-off. By summing all terms linear in SQ in the diagrammatic expansion 
of the free energy and taking the limit SQ — * one obtains the evolution equation for the free energy. Till now we 
have formally considered the total free energy, but for an inhomogeneous system it is useful to separate the bulk 



quantities, which refer to the homogeneous system without an external potential, from the surface terms due to the 
presence of the wall 



A Q = A? + AA Q 
p(r) = p b + Ap(r) 

F 2 Q (n,r 2 ) = Fg( Zl - z 2 , |si - s 2 1) + AF 2 Q (n,r 2 ) 
where the two point correlation function (ri, r 2 ) is the functional inverse of —C^i^i, r 2 ) i-e. 

J dv d F 2 Q (r 2 ,r 3 )C 2 Q (r 3 ,r 2 ) = -S(v 1 -v 2 ). (11) 

This differs from the known Ornstein-Zernike equation for inhomogeneous fluids because of the inclusion of the ideal 
gas terms. Note that, due to the cylindrical cut-off, F® b has lost the full rotational symmetry of the physical system 
for every Q G (0, 00). The separation in bulk and surface terms may be inserted in the diagrammatic expansion of the 
total Hclmholtz free energy. We begin with the derivation of the evolution equation for the bulk contribution to the 
free energy. Only the one loop diagrams containing 'bulk' F^ b bonds alone in the diagrammatic expansion contribute 
to the leading order in 8Q. The sum of these term can be carried out in closed form and reproduces the known 
random phase formal expression: 

0A^-^- SQ = -\ [p. [ ? f^ r log(l-#«, (12) 

2 J 2n J(Q,Q-SQ) ( 27r ) 

This equation is not yet suitable for taking the 5Q — ► limit because it contains the discontinuous function F® b (q, k) = 

— l/CH(qi, k). Expressing F^, in terms of the continuous function J^biQ' K ) = — K ) y i a Ecj s - © and (|11|) we 

obtain: 

iQ d - 2 J ^ log(l + Q)<f>(q, Q)) (13) 

where K&-\ is a geometric factor defined by 

^pr^- 1 (K-Q) = ^- 1 Q d - 2 . 

The initial condition must be imposed at Q — 00 and coincides with the mean field approximation for the Hclmholtz 
free energy of the fully interacting system. Unfortunately, equation (|13J) is not written in closed form because the 
evolution of the free energy depends on the knowledge of the two-particle correlation function ^"^(9, Q) at a generic 
value of cut-off Q. However, we can study how the two-particle function T^ h itself (or equivalently C^) is modified 
owing to the inclusion of fluctuations, by performing the same analysis leading to (|13fl . In this way we obtain an 
evolution equation for the pair correlation function, but this equation contains the three-particle and the four-particle 
correlation functions. The procedure may be iterated to higher orders thereby obtaining an infinite hierarchy of 
coupled differential equations. This hierarchy may be truncated at the level of the first equation by introducing a 
suitable anstanz for the two-particle correlation function and using the compressibility sum rule which relates the 
second density derivative of free energy to the zero wave-vector limit of the pair correlation function 

a 2 (-(3A?/V 
dp-. 




Q(k = 0)= v to '-■ (14) 



Note that this equation is exact for each Q-system, it provides an extremely useful link between thermodynamics and 
the long wavelength limit of correlation functions but it is not sufficient to close the hierarchy because it does not 
give information about the momentum dependence of T^iq, k). 

In the derivation of the bulk evolution equation, translational invariance allowed to perform the analytical resum- 
mation of the one loop diagrams leading to Eq. (|13f) without specifying the form of J 7 ^. Instead, in order to express 
in closed form the sum of the required diagrams in the expansion of the surface free energies, we need to specialize 
to a particular form of the inhomogeneous two point correlation function AF®(ri, r 2 ). Therefore, it is important to 
understand what type of fluctuations are described by AF^(ri, r 2 ). Let us consider the case of a wetting transition 



between two bulk phases, (3 and a, in contact with a wall. At a temperature larger than the wetting temperature and 
slightly off coexistence, the system will consist of a slab of the wetting phase, with density pp, an interfacial region and 
the bulk p a phase. If F£ b , which is the bulk correlation function of the a phase, is not equal to F^ b then AF®(ri, r%), 
which is defined as i*2( r ij r 2) — -^26 ( r i — r 2), describes not only the correlations in the interfacial region but also in 
the slab. In this case it is difficult to find a suitable ansatz for the form of AF^. However, in systems where the bulk 
correlation functions are the same for the two phases the problem simplifies and AF® just describes the interfacial 
region. This is precisely the case of the lattice gas model with nearest-neighbor interaction (i.e. of the Ising model) 
due to the symmetry pp = 1 — p a . If the inhomogeneous part of the two point correlation function describes only the 
interfacial region we can parametrize Ajft according to the standard factorized form [24| : 

A^( qu q 2 ,K) = g Q (qi)g Q (q 2 )F?(K) (15) 

where g®(z), the Fourier transform of g®(q), is usually identified as the first spatial derivative of the density profile 
and we have taken the wall perpendicular to z axis. Using this parametrization it is possible to perform the summation 
of the one loop diagrams which allows to obtain the evolution equation for the surface free energy: 

IQ i^P) = %V- 2 log(l + FiW(Q)) (16) 



where a® is: 



J 2tt l + 0(K,g)F 2 %,g) 



Note that the evolution of the surface free energy is coupled to the bulk evolution through the oft coefficient. 

In the general case of F£ b ^ F^ b the derivation of the evolution equation is more involved: By introducing further 
simplifying assumptions in the spirit of the local density approximation, it is still possible to obtain a decoupled 
evolution equation identical to Qlfijl with a slightly different coefficient oft |2fij : 

«2zm(«) = / ^ g Q {q)cp{n,q)g Q {~q). (18) 

Equation (|16f) has been obtained by starting from an expansion of the free energy at constant density profile: every 
Q-systems has the same p(r). It means that the external potential which stabilizes the density profile changes with 
Q. In the case of wetting or drying this procedure is clearly artificial: we would rather want to study how the density 
profile is modified by the inclusion of fluctuations at fixed external potential. In order to allow the change in p(r) 
when fluctuations are included, it is convenient to perform a Legendre transform of the total free energy jft . In this 
way, it is possible to obtain the evolution equations at fixed fugacity, i.e. by keeping fixed the quantity 

1 {z) = f3[p-U{z)], 

where U(z) is the microscopic interaction between the wall and the molecules of the fluid. We first define the (modified) 
grand free energy by 



-(3A Q + J d 3 r p Q {z) 1 {z) (19) 



Here 7(21) is fixed to the physical value, while the density profile changes with the cut-off wave vector Q. The 
Q-dependent density profile is related to w Q by: 



P Q (Z) = -ffy (20) 

As usual, for Q — » 00 oj® reduces to the mean field grand free energy and p®(z) approaches the density profile in 
mean field approximation. From the properties of the Legendre transform we obtain: 



9Q )~ \~dQ~), 



(21) 



As a consequence, the evolution equation for the grand free energy formally coincides with that of the Helmholtz free 
energy and the surface tension a® = Aw^ / S obeys the evolution equation: 



dQ 



-K d _ 1 Q d - 2 \og(l + FQ(Q)aQ(Q)). 



(22) 



Note that now a® is given again by Eq. (|1 T|) in terms of g®(q) and F^(k). Analogously to the bulk case, also surface 
quantities obey a hierarchy of coupled differential equations which can be closed at the free energy level ()22|l by a 
suitable ansatz for F®(k) and g®(z). Two surface sum rules will provide the link between the free energy and the 
correlation functions : 



dp 
d 2 aQ 
dp 2 



= T Q = dz 



P Q (z)-P?} = / dzApQ(z) 



J dz x dz 2 AT 2 Q (z 1 ,z 2 , k = 0) = f3F®{0) 



dz gQ(z) 



(23) 
(24) 



Equation 124|l can be identified as a "surface compressibility" sum rule while (|23|) constrains the form of the density 
profile. It is interesting to note that this two sum rules are sufficient to obtain a quite good description of the wetting 
phase transition [24| . 



B. Asymptotic Equations 

1. Bulk Equation 



In this section we show how the HRT evolution equations describe the asymptotic regime near phase transitions. 
The standard bulk HRT with spherical cut-off in the asymptotic critical region is known to give the same evolution 
equation obtained by momentum space RG approach Here we show that the introduction of a cylindrical cut-off 

does not modify this property of the bulk HRT equation. 

In the thermodynamic states close to the bulk critical point, equation (|13|l can be simplified. This regime is 
characterized by the growth of long-range fluctuations, that is by the divergence of ^F 2b (q = 0, n = 0) for Q — > i.e. 
in the final stages of the evolution. In order to study the critical region we have to extract the singular contribution 
to the free energy from the evolution equation (|13|) : Only a small neighborhood of the integration domain, close to 
q = and Q ~ contributes to the singularity. In this domain we can simplify the argument of the logarithm due to 
the divergence of the pair correlation function for q and n approaching zero. 




Q 



d-1 



dq 
2^ 



log 



Here qo is an arbitrary ultraviolet cut-off which is introduced to extract the singular contribution to the integral. Note 
that this equation does not contain the inter-atomic potential explicitly and acquires a universal form independent 
of the specific microscopic interaction. The next step is to choose a form for the bulk pair correlation function. As 
usual, we adopt an Ornstein-Zcrnike form for the direct correlation function 



c?& 0,<z) = 



-b(n 2 + q 2 ) + 



3_ 

dp 2 



V 



(25) 



where we imposed the sum rule Ijl4|l and we assumed that b remains finite also at the critical point. This particular 
ansatz for the pair correlation function implies that the critical exponent 77 is zero. We know that in three dimension 
7/ is not zero, but its value is small and it is zero to first order in the e-expansion pH so that, in three dimension, this 
is a good approximation. Substituting this parametrization in the evolution equation we obtain 




+ Q du 



Q 2 b(l + u 2 ) - — 




V 



We can eliminate the explicit dependence on the cut-off Q in the evolution equation by another change of variable 



i=-log(Q) 



Ht(z) = - 



A Q -A Q 



K d . 



where p^ c and Ab c are the density and the free energy at the critical point. Using these variables the asymptotic 
evolution equation becomes: 



dH t (z) 
dt 



= dH t {z) 



2-d 



zH' t {z) 



+00 



du 
2^ 



log 



H"(z) 



1 



H"(0) 



1 



where H'(z) and H"(z) are the first and the second derivative of H with respect to z and we set qo/Q 
Q — > limit. The integration can be easily performed and we obtain: 



dH t {z) 
dt 



= dH t (z) 



2-d 



1 



zH' t {z) + - [y/H»(z) + 1 - y/H' t '(0) + 1 



(26) 
co in the 

(27) 



This equation can be analytically investigated in d = 4 — e dimensions. The fixed point and the relevant eigenvalues 
provide the critical exponents to first order in e: 7 = 1 + e/6, (3 = 1/2 — e/6 and 6 = 3 + e which coincide with the 
known exact RG results j^. As expected, to lowest order in e, the choice of the cut-off symmetry does not modify 
the universal properties close to the bulk critical point within the HRT approach. 



2. Surface Equation 

The surface evolution equation we have obtained (|16|l depends on the form of the density profile through the quantity 
afi(tt) (|17I18|I if g®{z) is taken to represent the spatial derivative of p®{z). Here and in the following we assume that 
the shape of p®{z) is not affected by fluctuations. Fluctuations are assumed to shift rigidly the mean field density 
profile by a Q-dependent amount Iq identified as the distance of the interface from the wall: Ap®(z) = Apf(z — Iq). 
This assumption can be related to the non-homogeneity of the spectrum of fluctuations and is customary in most 
studies of the wetting phenomena [aUi|- From Eq.® we see that the fluctuations perpendicular to the interface, 
i.e. the fluctuations which distort the form of the profile, are only weakly divergent and as a first approximation we 
can neglect their effect. Note that this assumption is not justified near the bulk critical point, where the interplay 
between bulk and surface fluctuations does modify also the form of the density profile. Our analysis is therefore valid 
only away from the bulk critical point. The long wavelength behavior of the surface correlation function F^(k) is 
modeled by an Ornstein-Zernike form: 



-{Apf 



b,K z 



(28) 



where b s is a non-universal constant that tends to a finite limit at the wetting transition and we have used the exact 
sum rule (|2"Tjl : 



d 2 {-f3a c 



= (A p ) 2 F?(k = 0). 



(29) 



where Ap = J dzg^(z) = (p a — pp) is the difference between the bulk density of the two phases across the interface. 
Because we assume to be away from bulk critical point this quantity is regular at wetting transition. The assumptions 
we made on the density profile and on the surface correlation function are similar to those underlying the usual RG 
group treatment of the wetting transition 0, 0, IU El • Note that for wetting phenomena the critical exponent 77 is 
known to be zero P, 111 but the Ornstein-Zernike form for the surface correlation function is strictly accurate only in 
the interfacial region |18L fli^ . 

Away from the bulk critical point, the asymptotic form of the evolution equation for the surface free energy then 
becomes: 



dQ 



K, 



•g(-gqg)' 



KQ 2 



(30) 



where b' s — b s (Ap) 2 and we have disregarded non singular terms. Notice that ea. l|30|l does not depend on the specific 
shape of the density profile. The asymptotic evolution equation I j^Ojl is in fact fully equivalent to the non-perturbative 
"functional" RG approach jfll lldj. To prove this remarkable correspondence, we introduce the surface free energy 
a9 (r) by Legendre transform of a® (/i) : 



a Q (T)=a Q ( f j, Q )+fi Q T. 
Here T is the adsorption l)23|l and then [fe is such that 

da®(T) 



or 



(31) 



(32) 



The asymptotic evolution equation of aP(T) follows straightforwardly from the evolution equation of a® l)3()|) and 
from the properties of the Legendre transform 



d(-Pa9) _ K d ^, 2 ^J d 2 (-(3aQ) 



dQ 



-Q d -^og 



or 2 



b s Q' 



(33) 



By means of the change of variables: 



t = - 

i = r 



log(Q) 

K 



(d-3) . 

e 2 * 



we obtain the known RG equation 

dV t (l) 
dt 



v t (i) = % a® 



(d - l)V t (l) + i(3-d)^ + i log 



d 2 V t (l) 
DP 



(34) 



where Vt(l) is the effective interfacial potential in the RG approach. From this formal equivalence, we can identify 
the renormalizcd effective interaction Vt(l) as the surface Hclmholtz free energy a® and the height of the interface I 
as a measure of the absorption T. The initial (i.e. bare) form of Vt(l) is then given by the value of aP(T) at a suitable 
matching cut-off Q = Qq <c 1. The appropriate value should include the effects of short wavelength fluctuations (i.e. 
of fluctuations at Q > Qq) but a rough estimate is obtained starting from the mean field free energy. It is possible to 
show that, in the sharp-kink approximation |j| for the density profile of a continuous system, the initial form of our 
Vt(l) is given by (0 i.e. the same used in the effective capillary wave Hamiltonian. 



III. LATTICE GAS COMPLETE WETTING 



A. Mean field 



In this section we consider the wetting transition in a semi-infinite lattice gas model with nearest-neighbor interaction 
in contact with an unstructured wall. The interaction between the wall and the lattice gas is described by an external 
potential u,, where the subscript i > 1 labels layers of the lattice gas: we consider a planar geometry with an 
external potential translationally uniform in directions parallel to the wall. It is convenient to introduce an interlayer 
interactions which collects all inter-atomic nearest-neighbor interactions of the particles belonging to the z-th 
layer and the j-th layer, which are assumed to be nearest neighbors. Note that there is also a contribution va which 
describes the interaction of particles whitin the same layer. We shall discuss in the following only the particular model 
for which the external potential acts only on the first layer and is attractive 

Ui = uS^i 

and 

(v i,j nearest- neighbors layers 
4u i = j 
otherwise. 



with u, v < 0. In mean-field approximation, the free energy A m f of the system reads 

-(3A mf [3A R 1 



S 



(35) 



hi 



where pi is the density of the i layer, <f>ij = —(3vij, S is the area of the layers. The reference system is a hard core 
lattice gas. 



-PA R 



(36) 



It is well know that in mean-field approximation the bulk critical temperature is fcsT c /|t?| = 3/2. For T < T c there 
are two stable phases, a and /3, which coexist for /xq = — 3|u|. This value of the chemical potential at coexistence 
is not modified by fluctuations due to the particle-hole symmetry of the model, which also leads to a relationship 
between the bulk densities of the two coexisting phases: p a = 1 — pp. The mean field density profile is related to the 
external potential by thermodynamics: 



7» 



d ( (3 A 



mf 



d Pl V S 



= log 



1 - pi 



^2pj<t>n 



(37) 



where ji = (3{p — itj). This is an non linear equation for the density profile which can be conveniently written as: 

Pi = T^-T- (38) 



1 + e x 



where 



Xi = 7,: + },Pj4>ji- 



On general grounds, for nearest-neighbor interaction there are two possible scenarios depending on the value of 
u/v [26j named the strong substrate regime (u/v > 1) and the intermediate substrate regime (1/2 < u/v < 1). There 
is also another regime for < u/v < 1/2, the weak substrate regime, but its properties follow from those of the 
intermediate regime due to the symmetry of the Ising model: they describe the same phenomena of the intermediate 
regime but on "the other side" of the coexistence line pfij . i.e. in one case the bulk density is p a in the other case 
is pp. In FigQ] and Fig [21 we sketch the phase diagrams in the strong and intermediate regimes respectively. For 
m/w > 1 an infinite sequence of transitions occurs, corresponding to condensation of successive mono-layers. The 
critical temperatures of these layer transitions approach, for high-oder multi-layers, a well defined temperature T Rl 
the roughening temperature. At temperature below Tr isotherms for the adsorption T show an infinite sequence 
of sharp steps as p, — > fio- Between Tr and T c the steps are rounded but T still diverges approaching coexistence. 
Instead for u/v < 1 the layer transitions no longer extend to the T — axis but now meet the coexistence axis 
below a characteristic wetting temperature Tyy. Notice that the phase diagrams of the intermediate and the strong 
substrate regimes are similar for T > Tyy- 

Mean-field theory fails to describe the roughening transition and it incorrectly predicts Tr = T c |23|. Therefore, 
at mean-field level there is no smooth wetting transition but only a sequence of layer transitions, nevertheless it is 
interesting to note that in the intermediate regime for high temperatures the behavior of T(/i) at fixed T is very 
similar to that of a smooth wetting transition. We numerical solve the non linear mean field density profile equation 
H38J1 . with standard technique (2^, m the intermediate regime, u/v — 0.95, at different temperatures and we calculate 
the adsorption 

r = ^2(pi - Pb) 



where pt, is the bulk density obtained by 
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FIG. 1: Tipical surface phase diagram with representative gas-phase adsorption isotherms for strong substrate regime (u/v > 
is in units v and V is in units of lattice spacing. 
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FIG. 2: Tipical surface phase diagram with representative gas-phase adsorption isotherms for intermediate substrate regime 
l/2<u/u<l./iisin units v and T is in units of lattice spacing. 



At low temperature, as previously discussed, we see a sequence of layer transitions while at high temperature the 
layer transitions are still present but are shifted very near the coexistence line. As a consequence, the T(p) isotherm, 
shown in Fig [3 is smooth and looks very similar to what is expected in the wetting regime: The divergence close to 
the coexistence line is well represented by T(p) ~ A log (/i — /i co ) as in wetting. A fit of the data points gives a value 
A ~ 0.45 in a wide range of temperatures (see FigBJ. If we suppose that the Capillary Wave effective Hamiltonian 
describes the critical behavior of wetting phenomena, we can use our data on the critical amplitude A to estimate the 
wetting parameter ui: 

( 1 + l) = AA (39) 

where Ap is the difference between the bulk density of the two phases across the interface and we use I ~ T(Ap)~ 1 . 
Near coexistence Ap = 1 — 2pb due to the particle-hole symmetry of the system. The true correlation length ^| 
can be obtained from the Ornstcin-Zcrnike form of the bulk correlation we have adopted in this investigation: 

cosh^- 1 ) = I + i(£ 6 T 2 (40) 

(C b f - : P f B - ( P : )P ; (41) 
1 - 6pp b (l - p b ) 

where is the second-moment correlation length, analytically obtained by expanding to the second order the Fourier 
transform of the direct correlation function. By substituting into eq. (|39J) our numerical result for critical amplitude 
A, we obtain u> m f ~ as expected in mean field approximation. 



We can also evaluate the mean field surface susceptibility \ s defined as: 
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FIG. 3: Two absorption isotherm in the mean-field approximation for (5 = v/k^T = 0.72 and f3 = 0.9 (/3 C = 2/3) and u/u = 0.95. 
We see that for low temperature the layer transitions are evident, but for high temperature the isotherm is smooth, /j, is in 
units v and T is in units of lattice spacing. 
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FIG. 4: The logarithm plot of the adsorption isotherm in mean-field approximation for /3 = v/kbT = 0.74, (3 = 0.70, (3 = 0.68 
and u/v — 0.95. We see that in the asymptotic region the three adsorption isotherms are well approximated by straight lines. 
p, is in units of v and F is in units of lattice spacing. 



To obtain \ s in mean field, we differentiate the density profile equation i|37[) with respect to the chemical potential 



«9a 



- Pi) 



This is a tridiagonal set of linear equations in ^ that can be numerically solved. This route gives more accurate 
results than the direct numerical differentiation of the adsorption. Close to coexistence, the surface susceptibility 
diverges as Xs — A/A/i. 



B. HRT surface equation 



After having discussed the mean field approximation, which defines the initial condition of the HRT equation, we 
are ready to tackle the numerical solution of the HRT surface equation. In this first application of the HRT approach 
to inhomogeneous systems we are mainly interested in surface fluctuations in a temperature range not too close to 
the bulk critical point. In such a regime, the dependence of the surface evolution equation ill'L'l) on bulk quantities via 
the definition (|17[) can be neglected. Therefore, as first step, we have chosen to adopt a local density approximation 
for a Q (n) P|): 

a Q {n) = J g \gQ{q)\ 2 4>(q^). (42) 

In Appendix A we show how the evolution equation for the surface tension (I16fl , constructed for a continuum system, 
must be modified to be applied to a lattice system |2ll29|. The evolution equation for the surface tension of the lattice 
gas in three dimensions is given by 

d (-0a Q ) 1 

K q Q ' = 2^ 2(Q) l0g [1 + F 'W)«(Q)] (43) 

with Q G [—1 : 1] and Z?2(Q) is the two dimensional density of states (see Appendix A). 

According to the discussion in Section II-A, the factor g®(q) in equation (|42|l is related to the Fourier transform of 
the first derivative of the density profile. This definition, however, cannot be stricltly carried over to lattice models 
where the density profile is defined only on lattice sites. Instead, we parametrized g®(z) in terms of a rigidly shifted 
mean field density profile: 

gQ(z) = p' Q (x)=p' mf (z-l Q ) (44) 

where Iq is the (unknown) Q-dependent position of the interface and p m f(z) is an analytic interpolation of the mean 
field pi obtained by fitting the discrete density profile with a suitable continuum density profile: 

Pmf(z) = ^(a + b) + (a- b) erf (j^j ) . (45) 

This form does adequately represent the actual density profile by taking a equal to the bulk density and b to the 
contact density p\ . The width £j_ of the interface is only weakly dependent on the thermodynamic state and diverges 
logarithmically at wetting. As a first approximation we set £x independent of p. Instead, the position of the interface 
Iq does depend on the state and will be strongly renormalized by fluctuations. In fact, Iq is easily related to the 
adsorption Tq, via the parametrization (|44ll : 

Q (—.p,yQ\ r°° i 

dz {p m f{z-l Q )- p b ) ~—(p!- p b ) (46) 



dpp Jo yr ' ujy w/ r " J 2 

Note that, due to the parametrization (|44|l . the value of Iq drops out of the definition of oft(n) (|42|l which then 
becomes Q-independent (see Appendix B). We parametrized the surface correlation function by an OZ form: 

= P -4^ (47) 

JdzpUz) (l-A a«(/s)) 



in terms of the unkown parameter Xq. When this expression is substituted into Eq. 1)15(1 we get 

AFQ(z 1 ,z 2 ,K)=gQ(z 1 y 9 Q(z 2 ) (48) 

(1 - Aqot'(k)) 

where g Q (z) = p' Q {z)/\J dz' p' Q (z')\. The first term of (gEJ) can be interpreted as a normalized function which 
constrains the range of the correlation function to the interface, i.e. where p'q{z) does not vanish. The factor 
Pb(l — pb) guarantees the correct pb — > limit preserving the symmetry of the model. The chosen form of AF® 
(|48|l reproduces the momentum dependence of the known Random Phase Approximation |2^| in which oP{k) plays 
the role of an effective interaction potential on the interface. The yet unknown parameter Xq is determined via the 



surface compressibility sum rule l|24(l . valid for each Q-system. By use of the parametrization (|47ll in the sum rule 
B24JI we obtain the explicit relationship between Aq and the thermodynamics: 



d 2 {-Pa Q ) 
dW) 2 



dz p' Q (z) 



V = 0) = If-** (49) 
1 — AQa^(K = 0) 



Substituting the approximation of iv(ft) in the evolution equation (|43|l we obtain a closed non-linear partial differ- 
ential equation for the surface free energy <J^(^) at fixed temperature. The initial condition for this equation is given 
by the mean field theory. The boundary conditions are given at coexistence i.e. \x = \x co where we assume that the 
surface susceptibility diverges (wetting) and at \x — > — oo which corresponds to vanishing bulk density. The equation 
is written in quasi-linear form by a suitable change of variable shown in Appendix B. The evolution equation is then 
solved by a predictor-corrector algorithm 

In FigElwe show the logarithm of the susceptibility versus the logarithm of /i — /i co . Notice that the value of Xs 
at the last point of the grid depends on the mesh spacing as can be seen in FigEI F° r this reason this value is not 
considered in the following analysis. In the asymptotic region, close to the coexistence curve, log(xs) displays a linear 
behavior. We fit our data with a linear relation: \og(xs) ~ alog(/j, — fi co ) + b and we obtain a — — 1.040 ± 0.001. This 
value is compatible with a logarithmic divergence of the adsorption i. e. the inclusion of fluctuations does not modify 
divergence of the adsorption so, as predicted by RG approach, the critical exponent for complete wetting remain mean 
field like even when we introduce the fluctuations. The main effect of fluctuations in the asymptotic region is rather 
a renormalization of the amplitude A, see FigE] To estimate the value of the amplitude we fit the HRT results with 
log(xs) = log(A/A/i) = — log (a — ^i co ) + log A, see Fig0 We obtain A ~ 0.62 which is larger than the lattice mean 
field result A ~ 0.45. In Fig|8]we plot the surface susceptibility for two different temperatures. The fluctuations 
introduce a weak dependence of the critical amplitude on the temperature, which is absent in mean-field results. 




FIG. 5: The double logarithm plot of \ s isotherm for /3 = v/kbT = 0.72 and u/v = 0.95. We see that the effect of fluctuation 
is a renormalization of the critical adsorption amplitude. 

Using eq. (|39fl and our numerical result for A, we obtain ojhrt — 0.6 which is smaller than the RG estimation for the 
Ising model ui ~ 0.8 ^3 an d larger than the value u ~ 0.3 obtained from simulation results on critical exponents for 
critical wetting However our numerical estimate is obtained at higher temperature: T ~ 0.92T C . The deficiencies 
of the mean filed approximation prevents us to study a wide range of temperature and so it is also difficult to extract 
from our results how uj depends on temperature. However our data seem to imply an increasing of the wetting 
parameter for T — > T c (see. FigJSJ. 
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FIG. 6: The double Logarithmic plot of the Xa isotherm for j3 = v/ktT = 0.70 and u/v = 0.95 for two different numbers of 
mesh points. We see that the number of mesh points, i.e. the grid spacing, affects the value of Xs at the last point of the grid. 
The line is only a guide to eye 




FIG. 7: The double logarithm plot of Xs isotherm for f3 = v/kbT = 0.72 and u/v = 0.95 in the asymptotic region. The line is 
log(Xs) = _ l°g (a* _ Mco) + log A, with A = 0.62. 



IV. CONCLUSIONS 



In this work we have presented an extension of the HRT approach |2j to inhomogeneous systems and we applied the 
formalism to the wetting transition in the 3D Ising model in planar geometry. The approach presented here is, to our 
knowledge, the first microscopic theory able to describe capillary wave fluctuations at the level of the Renormalization 
Group approach. Its use is not restricted to the study of wetting but more generally to the study of fluctuations effects 
in non-homogeneous systems. Differently from the standard bulk HRT we introduce a cylindrical cut-off to take into 
account the anisotropic spectrum of the fluctuations in such a reduced symmetry. We obtain a hierarchy of evolution 
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FIG. 8: The double logarithm plot of \s isotherm for /3 = v/k b T = 0.70 and /3 = 0.74; u/v = 0.95. We see that the effects of 
temperature on the critical amplitude is small. The lines are only a guides to eye 

equations for the bulk properties which reproduces the known RG approach close to the bulk critical point together 
with an additional hierarchy of differential equations for the surface free energy and the interfacial correlations. This 
further hierarchy generally depends on bulk correlation functions. In order to close the hierarchy at the lowest level, 
i.e. keeping only the free energy equations, we analyzed a simple Ornstein-Zernike ansatz for both the bulk and the 
interfacial correlations. Within this approximation, the resulting bulk critical exponents turn out to be exact to first 
order in e = 4 — d and the surface free energy equation is equivalent to the known functional RG approach applied 
to the effective capillary wave Hamiltonian @, M, LjjJ] • To study the effects of fluctuations also in the non asymptotic 
region we investigated by the HRT approach the complete wetting transition in a lattice gas model with nearest- 
neighbor interaction (i.e. for the Ising model). In order to simplify the analysis we introduced a further decoupling 
approximation which disregards the effects of bulk fluctuations on surface quantities. The numerical results show 
an interesting renormalization of the critical adsorption amplitude which can be related to a renormalization of the 
wetting parameter. The HRT value of the wetting parameter is to ~ 0.6 which is smaller the the field-theoretical 
estimate for this model ^| but it is greater than simulation results for critical exponents of the critical wetting 
transition at lower temperature. This can probably related to the increase of wetting parameter for T — » T c . 
Notwithstanding the crudeness of some approximation we introduced, the results for critical and non-critical quantities 
compare quite favorably with available simulation data. Clearly there is room for improvement within the class of OZ 
closures of the HRT hierarchy. In particular it would be interesting to allow for a renormalization of the transverse 
correlation length £j_ which we kept constant. The extension of this approach to different systems, like off-lattice 
fluids or the Ising model in parallel plate geometry looks also promising. 
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APPENDIX A 



Consider a lattice system with nearest-neighbor interaction 




v i,j nearest-neighbors site 
otherwise. 



The Fourier transform of the dimensionless inter-particle potential is 

4>(q, k) = 4w/3[7 2 (k) + cos(g)] 

where w = —v > and 



1 d 

7d(«) = -T^cos(fci). 



d- 

i=l 

The Brillouin zone is defined by —n < fcj < n and we choose the cylindrical of the form [2J, |2_ 

W'^-\0 j 2 (k)>Q 
with Q G [—1,1]. The evolution equation for the surface tension of the lattice gas in three dimensions then becomes 

d(~8a Q ) 1 

dQ =2 D2{Q)l0g[1 + Fs{Q)am (50) 
where a(Q) is the value of aft(K) evaluated at 72 (ft) = Q and £>2(Q) is the two dimensional density of states: 

-■»(«))-; 

where A"(x) is complete elliptic integral of the first kind. 



APPENDIX B 



In this appendix we provide some detail on the algebraic manipulations on Eq. I|43f) necessary for implementing an 
efficient numerical algorithm. As a first step we give the explicit form for the effective surface potential obtained by 
substituting the parametrization discussed in the text for the function g®(z) 1(44(1 into the definition of o9(k) 142(1 . 
By use of the mean field density profile 1)45(1 g®(z) is just a Gaussian function and afi(K.) is independent of Iq: 



oft (k) ~ ^ b - Pl)2 dq e-i<> 2 &<t>(K, q) = l^—^l! f d q e~^2pw \{d - 1) 7 (*) + cos(g)] = 



2tt J_ v x ' 2tt 

= l3w( Pl - p h f [jd - 1)Ji 7 (k) + I 2 ] 

where 

h = I dqe-? q2 tt 



(51) 



/7T 
dq e~i q2 ^ 2± cos(g) 
-71 



Now we introduce the shorthand notation f r = pt, (1 — pi,) and 6q = —fta® . The surface susceptibility of the Q-system 
is easily expressed in terms of 6q : 

while the parameter Xq is given by eq. I(49|) . The evolution equation ((43(1 is written in quasi-linear form by introducing 
the new variable vq defined by: 

v Q = log (1 + P(Q)) 

where 

b^(p b - P i)- 2 fraoa(Q) b%f r r(Q) 



P(Q) = F„(Q)a(Q) 



6% (ao - a(Q)) + f r P 2 a(Q) 6£m(Q) + UP 2 a{Q) 



and 

vq = f r a{Q) + a 2 (Q)uQ. 

where a{Q) = a®(Q)(pb — pi)~ 2 ■ By differentiating twice the evolution equation (|43|l with respect to the chemical 
potential and substituting the definition of vq we get 

dV Q l D (0) d2vc > 



Inverting v — v(Vq) we find 



b" 



dQ 2 dp 2 

P 2 f r a(Q)[e v -1] 



Q [So - 5(Q)] [e» - 1] - / r a a(Q) 
Differentiating this relation with respect to Q we obtain 



d 



(3 2 f r a(Q)[e^ -1] 



9Q " dQ ^ [fio - d(Q)] [e-Q - f ] - / r a a(Q) 
where ''y '•" 



/3 2 / r 2 a(Q) 2 d t> o e^ - (i 2 f r aa [e^ - l] 
(m(g)[e«« -l] - f r r) 2 



(52) 



— *a — 9_vq_ nQw i n t roc [ uce the notation 



r(Q) = a &(Q) = (tj [(d - ifllQ + hh{d - 1)(1 + Q) + J 2 2 ] 
m(Q) = a - a(Q) = -(d- l)Ii(l - Q) 

7T 

which give f = aoa and m = —a. By solving the algebraic equation (|52|l for i)Q we obtain 

f r a\Q) 2f3 2 f 2 a 2 (Q)a w 

which is the evolution equation for vq . Finally, by a further change of variable, we introduce uq by 

f r a(Q) + &(Q) 2 u Q = log (1 + F s (Q)a(Q)) 



(53) 



(54) 



Switching from the chemical potential to the fugacity z 



_ p/80-Mco) 



we express the derivatives of vq in terms of uq: 



vq = f r a + 2aau Q + a uq 



v'A = a/ r ' + a 2 u Q = a/3 2 



9fr , 2 d2 fr 



dz ~^ dz 2 



a 2 (3 2 



du Q 2 d 2 u Q 



dz 



dz 2 



Substituting these expressions into the evolution equation for vq Ij53(l we obtain the quasi-linear equation satisfied by 
uq: 



du Q 
dQ 



Lx + L 2 +M— 

dz 1 



with 



Li = 


1 f " (e*« 
a 2 f r ct 2 


_ 1)V V « - / P 5 


L 2 = 


D 2 (Q)e- v i 


(m(Q)[e^-l] 


2(3 2 a i {Q)f 2 a 


M = 


D 2 {Q)e~^ 


(m(Q)[e*« - 1] 


2l3 2 a 2 (Q)f 2 a 



,dfr 

dz 



dz 2 



+ a 2 p 2 z 



duQ 
dz 
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